The Hall conductivity in unconventional charge density wave systems 
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Charge density waves with unconventional order parameters, for instance, with d-wave symmetry 
(DDW), may be relevant in the underdoped regime of high-T c cuprates or other quasi-one or two 
dimensional metals. A DDW state is characterized by two branches of low-lying electronic excita- 
tions. The resulting quantum mechanical current has an inter-branch component which leads to an 
additional mass term in the expression for the Hall conductivity. This extra mass term is paramet- 
rically enhanced near the "hot spots" of fermionic dispersion and is non-neglegible as is shown by 
numerical calculations of the Hall number in the DDW state. 
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Recently, the interest in charge density waves with un- 
conventional order parameters has increased 0, 0, El- 
la particular, it has been shown that a charge density 
wave with d-wave symmetry (DDW) represents a sta- 
ble state of the t — J model in the large-N limit in cer- 
tain doping and temperature regions 3]. It thus may 
be intimately related to the pseudogap phase of high- 
Tc superconductors^^!- The presence of a DDW state 
should also cause changes in transport coefficients^. 
Refs.0,0 discuss DDW-induced changes in the Hall ef- 
fect on the basis of the standard formulajg, which is ap- 
plicable to an usual metal. In the present paper we argue 
that a careful reconsideration of the Hall coefficient for 
the case of a DDW state results in an additional term to 
the usual expression. This term enhances the change in 
the Hall number due to onset of a DDW order parameter. 

Qualitatively, the appearance of the new term can be 
understood as follows. It is known from the band the- 
ory of metals that the quantum mechanical current op- 
erator consists of two parts, the intra-band derivative 
with respect to the wave vector and the term describing 
the inter band transitions. Usually, the inter band energy 
spacing is large enough to neglect the influence of the 
interband current term. 

In the DDW state the situation is different. The charge 
density wave with momentum Q couples electrons which 
differ in momentum by Q. Taking, for instance, a square 
lattice and Q = (tt, tt) a two-band picture is obtained in 
the reduced magnetic Brillouin zone. There exists regions 
around certain wave vectors K, the so-called "hot-spots" , 
where the quasi-particle energies of both bands are close 
to the Fermi level. The influence of these hot spots deter- 
mines the changes in the Hall conductivity, as shown in 
Ref.0. We find that the interband contribution to the 
current is particularly important in the vicinity of the 
hot spots, leading to significant changes in the theoreti- 
cal predictions. 

From a broader viewpoint, the necessity of inclusion 



of the interband current terms is known for the case of 
almost degenerate electron spectra and for the case 
of the electromagnetic response in nodal (d-wave) super- 
conductors. On a formal level, it can be illustrated as 
follows. Near the point of degeneracy, K, the spectrum 

can be represented as £k+K ~ \J + ^y- The inclu- 
sion of the external vector potential through the Peierls 
substitution, k — > k — eA, leads to a troublesome non- 
analyticity of the fermionic action on A, i.e. ^/|A| 2 . 
The recipe for the correct treatment of the electromag- 
netic response in such cases is well known It amounts 
to retaining the non-diagonal form of the Hamiltonian, 
which is analytic in A and contains the interband cur- 
rents, until the end of the calculation. 

The mean-field Hamiltonian in the DDW state is 



n 



kanaka- + iAka^ak+Q.o- + h.c. (1) 



Taking nearest and next-nearest neighbor hoppings t and 
t' into account and putting the lattice constant of the 
square lattice to unity, the electronic dispersion is £k = 
— 2£(cos fc^+cos k y )+4t' cos k x cos k y — n. In the following 
we also will use the abbreviations £± = (£k ± £k+Q)/2. 
The d-density wave the order parameter is of the form 
A k = Ao(cosfca; - cos k y ) = -A k+ Q. 

In terms of two-component fermion operator 

^ ko . = fflkCT' a k+Q cr) the Hamiltonian becomes 

w = E k , CT *LMk.with 



H k = 



6c «A k 



(2) 



It can be diagonalized by the unitary transformation U = 
expicr^k with 

9u = (1/2) arctan(A k /£_), 

where a 1 denote the Pauli matrices. We have UHU^ = 
diag{s\, £2) = h and the new quasiparticle energies are 
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e li2 = ± [£ + A k ] 



21 1/2 



(3) 
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The fermionic Green's function is given by 



G kty {iui) 



(4) 



and it is diagonalized by the same matrix U . We write 
G = U^gU with g = (uo — ft) -1 . The external vector 
potential is included into the Hamiltonian by the Peierls 
substitution — > -ffk-eA- 

In what follows we use the Kubo approach within the 
linear response theory. It expresses the d.c. conductivity 
tensor, a a p, in terms of the current-current correlation 
function, in the limit of static uniform external fields. 
|10| In case of point-like impurity scattering, allowing the 
neglectance of the vertex corrections to the corresponding 
diagrams, the standard derivation leads to the following 
formula 



, I dxdk dn(x) 
-e / ——^r — ^ — Tr 



(2tt) 3 dx 



V^G A (x)V^G R (x) , (5) 



where n(x) = (e x / T +l) -1 is the Fermi function, the sum- 
mation over the spin index has been performed, and the 
integration over k refers to the magnetic Brillouin zone 
in order to avoid double counting. Here G A ( R ' denote 
the advanced (retarded) Green's function defined on the 
real energy axis x. The group velocity, corresponding to 
the microscopic quantum-mechanical current, is 



(6) 



Using the above definitions the expression for the con- 
ductivity can be rewritten as Tr V a G A V p G R = 
Tr [v a g A vPg R ] with v a = WV a U and 

v a = d a h + [(U^d a U), h] = V a h. 

The last equation shows that in the basis which diagonal- 
izes the Hamiltonian, the current operator v a is defined 
by a "covariant" derivative, T> a , which differs from the 
usual derivative by the Christoffel symbol. Explicitly, the 
current operator is given by 



W l(2) 



dei 



(2) 



d Ai 



dk a 



(£ + Al )i/a dk a £_ 



(7) 



(8) 



The off-diagonal term in the above expression arises 
from the k-dependence of the unitary transformation U, 
and corresponds to the interband transition operator O. 
9] The "mass" operator in the new basis is d a d 13 HU . 
The explicit expression for it, 



V a V h 



8 2 h 
dk a dkp 



?7! a if 

Elk — £2k 



(9) 



is a smooth function in the whole Brillouin zone. 

The scattering processes are modelled by the imagi- 
nary part 7 of the poles of Green's functions, so that 




FIG. 1: Two diagrams contributing to the Hall conductivity 



9ii = (i^fiT V))- In the limit of a large scattering 
time, r = 7 -1 , the principal contribution to the conduc- 
tivity (jHJ is delivered by the combinations g A i(x)g R 1 (x) 
and g A 2(x)g R 2 (x), where the poles of the Green's func- 
tions differ only by the value for the damping. One eas- 
ily finds that this leading contribution to the conductivity 
contains only "intraband" velocity terms. In the limit of 
zero temperatures we have 



e 2 r 



dk 



H k ) 2 S(e lk ) + (1 - 2)] , (10) 



in accordance with previous findings Note that the 

"interband" current term V3 in the final expression for 
the conductivity is absent only in the d.c. limit, but is in 
general present in the optical conductivity tensor and also 
modifies the optical sum rule. The optical sum is defined 
by the new mass ©, averaged over the occupied states 
in the Brillouin zone, and should exhibit the deviations, 
- A 2 /E F in the DDW state. 

The next step is to evaluate the Hall conductivity ten- 
sor in the DDW state. The magnetic field can be included 
by considering the first-order change in the Green's func- 
tions due to the magnetic field in Eq. (JSJ , as discussed in 
|TT| . Writing B p = iA p x p and taking eventually the 
limit p — > 0, the change in the conductivity is described 
by the two diagrams shown in Fig. ^ 

The contribution from the left diagram in Fig. ^ as- 
sumes the form 



e^Tr 



V k °k V k+p/2 LT k+P V k+p/2 ( - T k 



(11) 



The second diagram is obtained from the above expres- 
sion by putting p — > — p and applying hcrmitian con- 
jugation corresponding to G A — > G R . The zeroth-order 
terms in p in Ea. Hll|) are odd in k and vanish upon the 
subsequent integration over k. The terms linear in p as- 
sume the form ^K^TyK^^ with the tensor K al3jr > 
given by 



g r v£g a 



r)V J 



r)C A 
qt>7 aLT k f/a 



f)V a 



n 



Further steps include the use of the property dG^/dk n = 
GkV^Gk, the application of the unitary transformation 
U with the corresponding change d a — > T> a , and an in- 
tegration by parts over k. Attention should be paid to 
the non-commutative property of the involved matrices. 
After some calculation K a ^ iri reduces to 



-W{g R v^){Vng A )v c 



g R v? 



[.</ 



(12) 
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momentum q x doping 



FIG. 2: The evolution of the Fermi surface with the opening 
of the DDW gap. 



Finally, we combine the expressions from the two dia- 
grams in Fig. ^ an d retain the principal contribution in 
the large-r limit. As a result we obtain for the Hall cur- 
rent j in the low-temperature limit 



r 



(13) 



^(eikXkWk 



/ d 2 e lk 
\dkpdk v 



£lk - £2k 



(14) 



with er-yn being the totally antisymmetric tensor. 

Eq. <|14|l is the central result of this paper. It shows 
that the Hall conductivity in the DDW state is defined 
by two inverse mass terms. The first term d 2 £i\Jdkpdk ri 
is the direct analog of the standard expression Q and is 
usually discussed 0,0- The second term ufk^k/ '( £ ik — 
£2k) is also present in but enters Eq. (|14(1 with an 
opposite sign. Let us discuss the relative importance of 
this term. 

First, this term contributes only in the anisotropic case 
and is unimportant particularly for an excitonic insulator 
[l2| , which is described by the Hamiltonian (|2J) with £k oc 
k 2 , £k+Q — k 2 and Ak = constant. In this case all 
three velocities, Vi,2,3 in Eq.©, are parallel to k. As 
a result, the second mass term in Ea. (|14fl . containing 
v i(2) x v 3: vanishes. 

Second, quite generally, the interband current is 
present for an electron in a periodic potential, so that 
the analog of (|14f) may occur in a multi-band metal as 
well. The main difference between this case and the dis- 
cussed DDW state lies in the relative importance of the 
second inverse mass term in l|14[l • The energy denomina- 
tor in it involves the interband splitting which is usually 
large in the multi-band case. The energies of two bands 
may become closer at the van Hove points in the Bril- 



FIG. 3: The zero-temperature results for the doping depen- 
dence of the conductivity a and the Hall conductivity <jh di- 
vided by e 2 r and e 3 r 2 , respectively. The contributions from 
the first and second inverse mass term in (1141 to an are shown 



(2) 



, H and a H , 



respectively. 



louin zone, however, the interband current v% vanishes 
there. These general arguments are inapplicable to the 
DDW situation as described below. 

For the anisotropic dispersion £k and DDW order pa- 
rameter Ak the second mass term in Eq. <|14|) is impor- 
tant. Indeed, the DDW-induced changes in the Hall con- 
ductivity are mostly determined by the vicinity of the 
"hot spots" in k— space where £k — £k+Q — 0. Expand- 
ing the spectrum around one of these spots we write 
i + ~ Viki, £- ~ V 2 k 2 , and A k ~ A hs + V d k x . Here 
jfei,2 = k x ±k y and |A hs | ~ \V d \ < |Vi| ~ \V 2 \. These ex- 
pressions lead to V3 ~ V 2 and a parametrically small en- 
ergy denominator £ik — £2k ~ |A| in l|14|) . Eq.l|14|l shows 
then an anomalously large contribution ~ V^V^/IAI in 
the hot spot's vicinity, 6k ~ IA0I/V2. Observing that 
T> a T)"h in © is finite near the hot spots, one expects 
that the second mass term enhances substantially the 
anomalous contribution from the first term in l)14|l . The 
resulting change in the Hall conductivity, 5o~ xyz , is esti- 
mated as 



e 3 r 2 27r- 1 I/ 2 y d si 5 n(yiA /(S ). 



(15) 



We see that 6a xyz is negative for the above form of the 
spectrum, thus enhancing the absolute value of the (neg- 
ative) a xyz . 

We emphasize that the correction Eq. (|15fl which is lin- 
ear in |Ao| explicitly contains the gap velocity Vd- In the 
case of an s-wave order parameter, Ak = Ao = constant, 
the velocity Vd = and the first nonvanishing correction 
to Sa xyz would be of order of A 2 , and thus much smaller. 

We have performed numerical calculations for a xx = a 
and o~ xyz — o~h using Eas. HlUfl and (|14|) and our Hamil- 
tonian Eq.Q. In rough agreement with the large-N 
limit of the t — J model [3] we modelled the gap by 
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FIG. 4: The temperature dependence of the Hall number 
n H = -(t 2 /cth- 

A = A(T)y(l+/i)e(l+/i), where fj, is the chemical po- 
tential, A(0) = 0.58, a BCS temperature dependence is 
assumed for A(T), and t is used as the energy unit. The 
onset of the gap at /i = — 1 corresponds to the critical 
doping 8 C ~ 0.145 at T = and to the critical tempera- 
ture T c ~ 0.064 at S = 0.075, using always t' = 0.3. Fig. 
121 shows Fermi lines of this model for three different dop- 
ings. The Fermi lines consist of arcs around the nodal 
direction and lines near the antinodal points. Lines for 
the same doping end at the boundary of the magnetic 
Brillouin zone at different points because of the presence 
of the gap. 

The conductivities at zero temperature were obtained 
as integrals over Fermi lines. We used several hundred 
points to parametrize the Fermi lines ensuring that sim- 
ilar grids were used for different lines to achieve a nu- 
merical cancellation of singular terms. The temperature 
dependent conductivities a(p,T) were calculated using 

f dn(x) f 1 
a(p,T) = J dx — — a(fi + x,0) — J dna(fi + x(n),0), 



with x(n) = Tln(n _1 — 1). The latter redefinition of 
the integration regularizes the calculation at low temper- 
atures. 

The conductivity a has a contribution linear in the 
order parameter coming from the vicinity of hot spots, 
5a xx ~ — e 2 T7T~ 1 \V2^hs/Vi\. It translates to a square 
root dip near the critical values 5 C and T c , as can be 
seen in the curve for a in Fig. 13 for the case of S c . As- 
suming that most of the scattering is due to impurities, 
t is qualitatively unchanged at T c . Consequently, the 
square root feature should be observable not only in ctr 
but also in a. Note, however, that this dip in an and a 
is determined by ^jr- and Ak at hot spots, respectively. 
As shown in Fig. y] the usual expression for an (first 
term in Ea.l|14|). denoted by o# ), exhibits only a very 
weak change at S c as a function of doping. In contrast 
to that, the new term (second term in Ea.ifTl|l. denoted 

(2) 

by a H ), shows a well-pronounced square-root behavior 
near S c and dominates the change in the total Hall con- 
ductivity <jh — + (7 i^- The temperature dependence 
of the conductivities is qualitatively similar to the doping 
one. 

Fig. 0] depicts the temperature dependence of the Hall 
number nn = —& 2 /&h- The curve denoted by is 
based on the usual expression, the curve n# on our com- 
plete expression including the extra mass term. The on- 
set of the DDW again causes an approximate square root 
decay below T c in both cases. From a quantitative point 
of view it is clear from this Figure that the conventional 
theory gives only roughly 2/3 of the decay so that the 
discovered new term cannot be neglected in quantiative 
calculations. 

In conclusion, we derived an expression for the Hall 
conductivity an in the CDW state including also the in- 
tcrband current contribution. As a result, there is an 
additional term to an which may be interpreted as a 
renormalization of the mass and which is especially im- 
portant for momentum-dependent CDW order parame- 
ters. It is shown numerically that the new term increases 
the anomalous contribution ~ \JT C — T to an by about 
a factor 2 in the case of the DDW. 
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